ISN 97.99 
DFTT 46/97 



Twist-4 effects and Q 2 dependence of diffractive DIS 



M. Bertini a , M. Genovesef] \ N.N. Nikolaev cd , A.V.Pronyaev 6 , 
: B.G. Zakharov d 

On 

•I— > 

^ ■ a INFN, Sezione di Torino, Via P.Giuria 1, 1-10125 Torino, Italy 

q \ b Institut des Sciences Nucleaires 

\ Universite Joseph Fourier-IN2P3-CNRS, 

53, avenue des Martyrs, F-38026 Grenoble Cedex, France 

^ . c IKP(Theorie), FZ Jiilich, 5170 Jiilich, Germany 

^ \ d L. D. Landau Institute for Theoretical Physics, GSP-1, 117940, 

O ' u l- Kosygina 2, Moscow V-334, Russia. 

t^. , Department of Physics, Virginia Polytechnic Institute and State University, 

0\ Blacksburg, VA 24061, USA 

a: 

Abstract 



X 

5-H 



In this letter we report the direct perturbative QCD evaluation of twist-4 effects in diffractive 
DIS. They are large and have a strong impact on the Q 2 dependence of diffractive structure 
functions at large (3. Based on the AGK rules, we comment on the possible contribution from 



diffractive higher twists to oc corrections to proton structure function at small x. These 
corrections to the longitudinal structure function Fl may be particularly large. 
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The direct calculation of higher twist corrections to the proton structure function (SF) 
remains one of challenges of perturbative QCD (pQCD) description of inclusive deep inelastic 
scattering (DIS) ([0, 0, for phenomenological estimates based on departures from the DGLAP 
evolution see ||). In a striking contrast to inclusive proton SF, twist-4 (T4) corrections to 
diffractive SF, F^ 3 \xjp, (3, Q 2 ), (i = L, T), are genuinely short distance dominated ||, [5| and 
thus calculable from the first principles of pQCD. Furthermore, they are strongly enhanced 
by the gluon SF of the proton squared f|, Fj D( ^(T4; xip, (3, Q 2 ) oc -^G 2 (xjp, Q 2 ). (Hereafter 
x,Q 2 ,(3 and xjp = § are standard diffractive DIS variables and G(x,Q 2 ) = xg(x,Q 2 )). For 
instance, for (3 > 0.9, the twist-4 longitudinal diffractive SF was found to be comparable to, 
and even exceeding, the leading twist-2 (T2) transverse SF over a broad range of Q 2 , making 
questionable || [7], §, |J applications of the standard DGLAP evolution to diffractive SF. 

In this letter we extend our early analysis of twist-4 longitudinal diffractive SF |J to 
similarly enhanced twist-4 contribution to the large-/3 transverse diffractive SF |. It enters 
F^ 3 \xjp, (3, Q 2 ) with negative sign and, as the higher twist dies out at large Q 2 , the transverse 
SF (xjp, (3, Q 2 ) will rise with increasing Q 2 . We present for the first time an analysis 
of the (3 and Q 2 dependence of compensations of the longitudinal and transverse twist-4 
SF's in = F®^ + F^ 3 \ which completes the evaluation of driving twist-4 terms in 

F^ 3 \ We compare our results with the HI experimental data on F^ 3 ^ \T1\. Based on a 
connection between diffractive DIS and unitarity corrections to the proton SF suggested by 
the Abramovski-Gribov-Kancheli (AGK) rules |T2|, 13], we evaluate the unitarity-driven 
diffractive twist-4 corrections to the small- £ proton SF, which is potentially very large in the 
proton longitudinal SF Fl- 

The underlying pQCD subprocess for diffractive DIS at large (3 is the excitation of the 

qq Fock states of the photon 7*p — > Xp', where X = qq ||. The relevant pQCD diagrams 

are shown in Fig. 1. In what follows, z and (1 — z) are the fractions of the (light-cone) 

momentum of the photon carried by the quark and antiquark, respectively, k is the relative 

transverse momentum in the qq pair, M 2 = (m 2 + k 2 )/z(l — z) is the invariant mass of the 

diffractive system, p± is the (p,p') momentum transfer and t = —p± 2 . The contribution of 
2 The preliminary results from this study have been reported at DIS'97 ||. Related results were presented 
at DIS'97 by J.Bartels and M.Wiisthoff @. 
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excitation of one open flavour of electric charge e/ (in units of the electron charge) and mass 
in f to transverse (T) and longitudinal (L) diffractive SF's for forward diffraction dissociation 
is calculable in terms of the helicity amplitudes $1 and $2 introduced in |J: 



F° i4) (t = 0,XTP,f3,Q 2 ) 



Aire 2 (3 r dk 2 (k 2 + m 2 f ) 



*UQ 2 ) 



k 2 +m 
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where J = y 1 — 4 f M2 is the Jacobian peak factor []. To the Leading Log^ approximation 
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Here f(x, k 2 



{k 2 + m 2 ) 2 {k 2 + m 2 ) 3 
dG(x,k 2 )/d\og(K 2 ) is the unintegrated gluon SF of the proton, e 2 = z(l — 



z)Q 2 + mj, a = e 2 + k 2 + 



K 



2kK, 



r 2 , . 2 _ k 2 + mj 
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is the pQCD scale of the Leading LogQ 2 approximation (J7], [RJ), and we also have shown the 
Leading LogQ 2 approximations for $1, $ 2 . Q 2 has already been used as such in flip, (0) as the 
argument of the strong coupling as- The above outlined formalism has become a standard 

7], O, I5L ITJi]. Extra contributions, oc f(xjp, Q 2 ), which take 



description of DIS at large (3 



over in $1, <E>2 at /? <S 1, were found in fL6[ , but they can be neglected at (3 ~ 1. 

3 We follow the convention in which diffractive structure functions are dimensionless hence 

vtotipp) = 40mfe in the denominator of the r.h.s of Eqs. (]l|), (||). Our definition (|^) of F®^ differs from 
O, by the factor xip, so that our does not blow up ~ -p- when xjp — > 0. 



The advantage of forward diffraction is that the SF F^ 4 \t = 0, xjp, f3, Q 2 ) is calculable 
directly in terms of the conventional gluon SF. The first measurements of the t-dependent 
F^\t, xjp, (3, Q 2 ) have been reported recently by the ZEUS collaboration [n|, but for large (5 



of our interest in this study only the t- integrated diffractive SF's are available ||11|| . Assuming 
the usual da/dt oc exp(— Bdp 2 ^), one finds 

F^\ X v,f3,Q 2 )= fdt a ^F^\t, X v,(3,Q^ (5) 
J lo7T Lbirn^p) 

Using Eqs. (Q),(^), one readily finds that the longitudinal SF F®^ is dominated by the 
contribution from large k 2 ~ \M 2 , which is further enhanced because of scaling violations 
in the gluon structure function of the proton. Consequently it is perturbatively calculable 
||, ||, |6| and at large Q 2 one finds the higher twist []. 

F? (3) (T4; t = 0,x r ,fi, Q 2 ) = AP3 H ~ W? ■ Q 2 ) , (6) 



Q 2 5 3 ip 12 



where 



^(xip,/3,Q 2 )= /^f 1 -^) 2 [a s (Q 2 )G(xjp,Q 2 )f , (7) 



1-/3 



and for quick estimates to a logarithmic accuracy, J r g(xjp, j3, Q 2 ) ~ [as(\Q 2 )G(xjp, \Q 2 )} 2 - 

For the separation of the transverse SF F^^ into the twist-2 and twist-4 components, we 
notice that 

k 2 + m 2 (3 k 2 + m 2 

M 2 = T^P ' Q 2 ' (8) 

which defines the decomposition 



^>(T4; . = 0, ft Q ) = - 6(1 _^ Ba , / ' <4(<? )*i ■ (10) 



4 We emphasize that in this paper we focus on (3 ~ 1. As it has been shown in ]l7[ , at /3 <C 1 , dominated 
by the ggg excitation, the longitudinal diffractive SF is twist-2 and R = ul/&t ~ 0.2, in close similarity to 
inclusive DIS M. 
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Notice, that the large-A; 2 behaviour of the integrand of the the twist-4 transverse SF F^^\TA) 
is identical to that of i^jP , it is pQCD calculable, and for large Q 2 , well above the charm 
threshold, we find 



F^ (3) (T4;t = 0,Xv>,P,Q 2 ) 



In contrast to that, the twist- 2 transverse SF is dominated by k ~ mj : 



F^ {4) (T2;t = 0,X!P,/3,Q 2 ) oc 



dk 2 k 2 



m 2 + k 2 1 ~ 



(12) 



(k 2 + m 2 f ) 3 

(we suppress in QT2] ) terms containing extra factors [m 2 /(&; 2 + m 2 )]™). Still, for heavy flavours 
and/or sufficiently small (1 — (3), the pQCD scale Q 2 is large, one is in the legitimate pQCD 
domain and to a logarithmic accuracy, 



F° {3 \T2-x F ,P,Q 2 ) 



/ 3(l-/3) 2 (3 + 4/? + 8/3 2 ) e) 



Qm 2 B d {T2-f3) 



12 



as{Q 2 )G 2 (x w ,Q 2 



m 



f 



.(13) 



^(4) 



Both terms, oc $ 2 and oc $| i n Eq. ©, give comparable contributions to the twist- 2 F^' 
Because of the finite phase space of the k 2 integration, k 2 < \M 2 — m 2 , there emerge higher- 
twist looking small corrections oc ( jfi 1 , n> 1, to the above defined twist-2 SF. These 
corrections are logarithmically weaker than the genuine twist-4 transverse SF ( JTT| ) and for 
this reason, we do not separate them. 

Evaluating the twist-2 transverse SF at not so large f3, one needs a model for the small-Q 2 
behaviour of the unintegrated gluon structure function f(x,Q 2 ). Because of cancelation of 
soft gluon radiation from different quarks in the colour singlet nucleon, f(x, Q 2 ) vanishes at 
Q 2 — > 0. The Born approximation suggests the soft Q 2 behaviour [ l6|, pfl 

dG B (Q 2 ) AC B a s (Q 2 ) 



[l-G 2 (Q,-Q)}, 



(14) 



dlogQ 2 7T (Q 2 + fi 2 G ) 2> 

where G2{Q, —Q) is the two-quark form factor of the nucleus, G(0, 0) = 1. For harder gluons 
one must interpolate between the form ( |i"4"|) and any convenient parameterization for G(x, Q 2 ), 
for instance, the GRV NLO parameterization [p2]| : 



G(x, Q 2 ) = G B {Q 2 



' Ql 
Ql + Q 2 



N 



G 



GRV 



>,Q 2 



' Ql ' 
,Ql + Q 2 , 



N 



(15) 



For a similar soft Q 2 parameterizations see Ref. [|T^]. The GRV formulas hold only for 
Q 2 > Ql = 0.4 GeV 2 , we take Ggrv{x,Q 2 < Q 2 ) = Ggrv{x,QD- We choose the parameters 
C B = 1.5, fi G = 0.15GeV, Q 2 = 3GeV 2 , N = 1 so as to reproduce the colour dipole cross 
section ||23|| . This way, our results for diffractive DIS are very close to those found in the 



colour dipole model |fL7| , 0]. For the quark masses we take m c = 1.5 GeV, m s = 0.3 GeV and 
mu,d — 150 MeV as in |17 |. For light flavours, Eq. (13) must not be taken literary, in this 
case the pQCD scale Q 2 is set not by m 2 but rather by a range of rapid variation of the gluon 
structure function, Q 2 ~ 0.5 GeV 2 . 

In order to calculate we need to know the diffraction slope B d . It has been anticipated 
some time ago ||, |5|, |TJ, |T7| that the typical scale for the diffraction slope B d (/3) is set by 
-B31P ~ 6 GeV -2 , as seen in hadronic diffraction in the so-called triple pomeron region, which 
in DIS corresponds to (3 <C 1. This prediction has been confirmed by the first data from the 
ZEUS LPS: B d = 7.1 ± l.lt% GeV -2 in diffractive DIS for 5 GeV -2 < Q 2 < 20 GeV -2 g§ 
and B d = 7.7±0.9±1.0 GeV -2 in real photoproduction |p5 |. For the numerical calculation of 



the twist-2 transverse SF we will use the results of Ref. || [26R , where it has been calculated 
that Bd(T2) w B31P for (3 — > 0, it rises by ~ 50% reaching a maximum at (3 ~ 0.5, then drops 
back to -B31P at (3 ~ 0.9, has a sharp minimum at (3 — > 1 and ends up at Bd((3 = 1) ~ -B31P 
in the exclusive limit of vector meson production. For the short-distance dominated twist-4 
SF's only the proton size contributes to the diffraction slope Bd and, consequently, we will 
use B d (T4;{3)^ B 3JP 

For comparing the large-/? behaviour of the twist-4 and twist-2 contributions to the trans- 
verse SF at large /3, it is sufficient to use the logarithmic approximation, in which 



F° (3 \T4 ] t = 0,x F ,f3,Q 2 



F° (3) (T2;t = 0,x F ,(3,Q 2 ) 



48m 2 



B d (T2;P) 



(3 + 4(3 + 8(3 2 )Q 2 £? 3 ip 







a s (\Q 2 )G(x F ,lQ 2 
a s (Q 2 )G(x W) Q 2 



(1-/3)' 



(16) 



The rise of this ratio oc when (3 — * 1 is to a large extent compensated by the growth of 
G 2 (xjp, Q 2 ) because of the rising pQCD scale Q 2 oc yz^ and, as a minor effect for (3 > 0.9, by 
the decrease of B d (T2; (3) with (3 — > 1 as it was found in [26 |. 
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If the longitudinal polarization of the virtual photon equals e^, then the twist-4 T and L 
SF's enter in the combination F2 tL = Ft + ^lFl- To the leading LogQ 2 approximation, one 
readily finds 



■2 



F2 3) (T4; xn>, P, Q 2 ) = T^ns—^K 1 + 2e L )(l - 2(3) 2 - l]^(a*, (3, Q 2 ) , (17) 



and the overall twist-4 contribution is dominated by the the positive valued for (3 beyond 
the crossover points, (3 > (3 + or (3 < and by the negative valued F^ A \TA) in between, 
where 

*= = 5 1 1 ± tt^t) (18) 

and (3 + = 0.79 for = 1. Because the subleading contributions FjP become larger as (3 
decreases, cf. Eqs. (|[) and @, the crossover is shifted to (3 + ~ 0.65. 

In Fig. 2 we show our results for F®^ and its decomposition into the transverse and 
longitudinal (assuming = 1), and the twist-2 and twist-4 components for xjp = 1.33 10~ 3 . 
Because of the inequality for the excited mass, M 2 > 4m 2 open charm excitation contributes 
only at Q 2 > Am 2 -^ and the charm threshold effects are clearly visible in Fig. 2. As the 
suppression oc ^2 of the open charm contribution to transverse twist-2 SF is to a large extent 
compensated by the large pQCD scale Q 2 , the charm abundance is substantial even in F^ 3 \ 
It is still larger in the longitudinal SF, in which the flavour dependence only enters via the 
lower integration limit in Eq. (^). The fractional contribution of open charm to F^ (3) (T4) 
is similar to that in F^ 3 \ but the charm threshold effects are smoother because $J oc k 2 
which suppresses the near-threshold cross section. Apart from the charm threshold effects, 
the twist-2 transverse SF is practically flat vs. Q 2 . Whereas for (3 = 0.65 the longitudinal 
and transverse twist-4 contributions cancel each other, at — 0.9 the twist-4 contribution is 
very large. Notice how scaling violations slow down the rate of decrease of twist-4 structure 
functions at large Q 2 . We find good agreement between our results and the recent experi- 
mental data on F®^ at xip = 1.33 10~ 3 from the HI collaboration [jTT]|. In Fig. 3 we show 



our results for the x^p dependence of F^(xjp, (3, Q 2 ) in comparison with the HI experimental 
data. Our formalism describes excitation of the continuum states, and the real comparison 
between theory and experiment is justified only for masses above the IS and 2S vector meson 
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excitation, M 2 > 3GeV 2 , which for (3 = 0.9 implies Q 2 = M 2 ^ > 20-30 GeV 2 . In the 
spirit of the exclusive/inclusive duality ||, our results for the continuum in the vector meson 
excitation region must be applicable to the experimental data in which the vector meson 
peaks have been smeared out. With these reservations, our calculations reproduce rather well 
the experimentally observed Q 2 and xjp dependence of F^ 3 \xjp,f3, Q 2 ) over the whole range 
of Q 2 and xip. Notice the steeper xjp dependence at smaller xjp, which is a characteristic 



prediction of pQCD models of diffractive DIS |L7| [18], |19|1 . The accuracy of the experimental 
data on diffractive DIS is not yet sufficient to test different models for gluon density in the 
proton. Theoretical evaluations of the short- distance dominated and F^ 3 \T4) must 

be good to ~ 10-20% accuracy (due to uncertainties in G(x, k 2 ) in the perturbative region), 
whereas for the leading twist F®^ 1 (where one has a large penetration into the low-fc 2 region) 



a conservative estimate for the accuracy is ~ 20-30% [17] . 

It is interesting to comment on the possible impact of diffractive higher twist on the 
proton SF at small x. The experimentally observed small- x growth of the proton structure 
function can not go indefinitely without running into conflict with the unitarity. In general, 
one can decompose the experimentally measured proton SF into the Born, alias the DGLAP, 



component and the unitarity correction [13], [14]] : 



F 2 (x, Q 2 ) = F 2 (DGLAP; x, Q 2 ) + A u F 2 (x, Q 2 ) (19) 

The standard DGLAP evolution, which is linear in parton densities, is strictly applicable only 
to the non-unitarized SF F2(DGLAP;x,Q 2 ). The unitarization of the proton SF is one of 
the open non-perturbative problems. The AGK rules suggest (i=T,L,2) [|1^, [L3[ 

A c/ F l (x,g 2 ) = -|_^^ D(3) (x IP ,AQ 2 ), (20) 

where x max =0.05-0.1 is the threshold for diffractive DIS. We are interested in the twist-4 
component of the unitarity correction coming from the diffractive twist-4 component of qq 
excitation. Because the integral (|20|) is dominated by the contribution from (3 ~ 0.5, one can 
put (3 = 0.5 in the scaling violations factor (^) and take it out of the integral (|20|) . Then, 
after little algebra, one finds 



S 



e 2 

AuFiiT^x^Q 2 ) = / Jg(x,0.5,Q 2 ) 



3, for i=T 
-4, for i=L (21) 
-1, for i=2 

which does not depend on x max as soon as x <C x max . It is convenient to present the diffractive 

Q 2 

twist-4 effects as the correction factor (1 — ^r) to the leading-twist structure function. Our 
results for the so defined Q\ 2 are shown in Fig. 4. We find that Q\ 2 rises with decreasing x 
and increasing Q 2 , but it is relatively small because of the strong cancelation of the longitu- 
dinal and transverse twist-4 contributions. It is interesting to notice that diffractive twist-4 
correction to the longitudinal structure function of the proton is much larger, 

QIl = |<24, 2 ~ 20Q 2 i2 , (22) 



where we use R = <?l/o~t ~ 0.2 PU| (for an indirect experimental estimates see 

Of course, the above estimated diffractive higher twist from diffractive DIS at large-/? is only 
one of the many possible sources of twist-4 corrections to the proton structure function at 
small x. We focused on it for the reasons that it is parameter-free calculable in pQCD, is 
enhanced by the gluon structure function of the proton squared and receives a substantial 
contribution from the heavy flavour excitation. At the moment, one can not exclude that a 
comparable contribution to Q\ i will come from the small-/? diffraction. 



Summary and Conclusions. 

We have presented the calculation of the Born approximation for the higher twist-4 correc- 
tions to the transverse diffractive structure function of the proton. It is similar to the longitu- 
dinal diffractive structure function which to the Born approximation is entirely twist-4. Both 
twist-4 diffractive structure functions are short-distance dominated and are perturbatively 
calculable. Their contribution to the diffractive structure function is large and dominates 
F®^ at (3 > 0.9. Evidently, this twist-4 component of must be removed before ven- 

turing the DGLAP evolution analysis of the experimentally measured diffractive structure 
function. Because the pQCD calculation of the twist-4 corrections to F®^ is parameter free, 
such a correction of the diffractive structure function must not cause any ambiguities. 
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Taking for guidance the AGK rules, we also evaluated the impact of the diffractive twist-4 
to the twist-4 terms in the proton structure function. The estimate (|22| ) serves as an warning 
that diffractive twist-4 corrections to small-x longitudinal structure functions of the proton 
can be surprisingly large. 
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Figure captions 

Figure 1 : One of the 16 Feynman diagrams for diffractive cross section. 

Figure 2 : The Q 2 dependence of the diffractive structure function and its decomposition into the 
transverse (T) and longitudinal (L) components and into the twist-4 (T4) and twist- 
2 (T2) components. The open charm threshold in (xip, /3, Q 2 ) is seen from the 
comparison of solid curves for excitation of 4 flavours (u + d + s + c) and dotted curves 
for excitation of light flavours (u + d + s) . The predictions are compared with the HI 
experimental data on F^ 3 \xjp, (3, Q 2 ) at xip = 1.33 10~ 3 [|ll|] for (3 = 0.65 (box a) and 
(3 = 0.9 (box b). 

Figure 3 : Our predictions for the xjp dependence of diffractive structure function for /3 = 0.65 
(the two top raws) and f3 = 0.9 (the two bottom raws) in comparison with the HI 



experimental data |]TTJ for several values of Q . 



Figure 4 : The scale Q 2 2 for the AGK evaluation of the diffractive higher twist contribution to 
the proton structure function at small x. 
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Figure 1: One of the 16 Feynman diagrams for diffractive cross section. 
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p=0.65 



(a) (3=0.9 



(b) 




Figure 2: The Q 2 dependence of the diffract ive structure function and its decomposition into 
the transverse (T) and longitudinal (T) components and into the twist-4 (T4) and twist-2 
(T2) components. The open charm threshold in F^ 3 \xsp, j3, Q 2 ) is seen from the comparison 
of solid curves for excitation of 4 flavors (u + d + s + c) and dotted curves for excitation of 
light flavors (u + d + s). The predictions are compared with the HI experimental data on 
F° {3) (xjp,p,Q 2 ) at xyp = 1.33 1(T 3 |1 for p = 0.65 (box a) and (3 = 0.9 (box b). 
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Figure 3: Our predictions for the xip dependence of diffractive structure function for (3 = 0.65 
(the two top raws) and f3 = 0.9 (the two bottom raws) in comparison with the HI experimental 
data [II for several values of Q 2 . 
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